Infinitesimal local operations and differential conditions for 

entanglement monotones 

Ognyan Oreshko\Q 

Department of Physics, University of Southern California, Los Angeles, CA 90089 

Todd A. Bruifl 

Communication Sciences Institute, University 
of Southern California, Los Angeles, CA 90089 
(Dated: 1st February 2008) 

Abstract 

Much of the theory of entanglement concerns the transformations that are possible to a state 
under local operations with classical communication (LOCC); however, this set of operations is 
complicated and difficult to describe mathematically. An idea which has proven very useful is that 
of the entanglement monotone: a function of the state which is invariant under local unitary trans- 
formations and always decreases (or increases) on average after any local operation. In this paper 
we look on LOCC as the set of operations generated by infinitesimal local operations, operations 
which can be performed locally and which leave the state little changed. We show that a necessary 
and sufficient condition for a function of the state to be an entanglement monotone under local 
operations that do not involve information loss is that the function be a monotone under infinites- 
imal local operations. We then derive necessary and sufficient differential conditions for a function 
of the state to be an entanglement monotone. We first derive two conditions for local operations 
without information loss, and then show that they can be extended to more general operations by 
adding the requirement of convexity. We then demonstrate that a number of known entanglement 
monotones satisfy these differential criteria. Finally, as an application, we use the differential con- 
ditions to construct a new polynomial entanglement monotone for three-qubit pure states. It is 
our hope that this approach will avoid some of the difficulties in the theory of multipartite and 
mixed-state entanglement. 
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I. INTRODUCTION 



One of the most impressive achievements of quantum information theory is the theory of 
entanglement. The theory of entanglement concerns the transformations that are possible 
to a state under local operations with classical communication (LOCC). The paradigmatic 
experiment is a quantum system comprising several subsystems, each in a separate labora- 
tory under control of a different experimenter: Alice, Bob, Cara, etc. Each experimenter can 
perform any physically allowed operation on his or her subsystem — unitary transformations, 
generalized measurements, indeed any trace-preserving completely positive operation-and 
communicate their results to each other without restriction. They are not, however, allowed 
to bring their subsystems together and manipulate them jointly. An LOCC protocol consists 
of any number of local operations, interspersed with any amount of classical communication; 
the choice of operations at later times may depend on the outcomes of measurements at any 
earlier time. 

The results of Bennett et al. [l|, 0, 0] and Nielsen 01 , among many others [1, H, 0, IE 0| > 

have given us a nearly complete theory of entanglement for bipartite systems in pure states. 
Unfortunately, great difficulties have been encountered in trying to extend these results both 
to mixed states and to states with more than two subsystems {multipartite systems). The 
reasons for this are many; but one reason is that the set LOCC is complicated and difficult 
to describe mathematically jl^ . 

One mathematical tool which has proven very useful is that of the entanglement mono- 
tone: a function of the state which is invariant under local unitary transformations and 
always decreases (or increases) on average after any local operation. These functions were 
described by Vidal and large classes of them have been enumerated since then. 

We will consider those protocols in LOCC that preserve pure states as the set of opera- 
tions generated by infinitesimal local operations: operations which can be performed locally 
and which leave the state little changed including infinitesimal local unitaries and weak 
generalized measurements. In Bennett et al. it was shown that infinitesimal local oper- 
ations can be used to perform any local operation with the additional use of local ancillary 
systems-extra systems residing in the local laboratories, which can be coupled to the sub- 
systems for a time and later discarded. Recently we have shown that any local generalized 
measurement can be implemented as a sequence of weak measurements without the use of 
ancillas This implies that a necessary and sufficient condition for a function of the state 
to be a monotone under local operations that preserve pure states is the function to be a 
monotone under infinitesimal local operations. 

In this paper we derive differential conditions for a function of the state to be an entan- 
glement monotone by considering the change of the function on average under infinitesimal 
local operations up to the lowest order in the infinitesimal parameter. We thus obtain con- 
ditions that involve at most second derivatives of the function. We then prove that these 
conditions are both necessary and sufficient. We show that the conditions are satisfied by a 
number of known entanglement monotones and we use them to construct a new polynomial 
entanglement monotone for three-qubit pure states. 

It is our hope that this approach will provide a new window with which to study LOCC, 
and perhaps avoid some of the difficulties in the theory of multipartite and mixed-state 
entanglement. By looking only at the differential behavior of entanglement monotones, we 
avoid concerns about the global structure of LOCC. 

In section II, we define the basic concepts of this paper: LOCC operations, entanglement 
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monotones, and infinitesimal operations. In section III, we show how all local operations 
that preserve pure states can be generated by a sequence of infinitesimal local operations. In 
section IV, we derive differential conditions for a function of the state to be an entanglement 
monotone. There are two such conditions for pure-state entanglement monotones: the first 
guarantees invariance under local unitary transformations (LU invariance), and involves only 
the first derivatives of the function, while the second guarantees monotonicity under local 
measurements, and involves second derivatives. For mixed-state entanglement monotones 
we add a further condition, convexity, which ensures that a function remains monotonic 
under operations that lose information (and can therefore transform pure states to mixed 
states). In section V, we look at some known monotones-the norm of the state, the local 
purity, and the entropy of entanglement-and show that they obey the differential criteria. 
In section VI, we use the differential conditions to construct a new polynomial entanglement 
monotone for three-qubit pure states which depends on the invariant identified by Kempe 
[l3l |. Finally, in section VII we conclude. In the appendix, we show that higher derivatives 
of the function are not needed to prove monotonicity. 



II. BASIC DEFINITIONS 
A. LOCC 

An operation (or protocol) in LOCC consists of a sequence of local operations with 
classical communication between them. Initially, we will consider only those local operations 
that preserve pure states: unitaries, in which the state is transformed 

p UpU\ U^U = UU^ = J, (1) 

and generalized measurements, in which the state randomly changes 

p^p,=M,pM]/p„ Y.M]M,=i, (2) 

3 

with probability pj = Tr |MjMjp|, where the index j labels the possible outcomes of the 

measurement. Note that we can think of a unitary as being a special case of a generalized 
measurement with only one possible outcome. One can think of this class of operations as 
being limited to those which do not discard information. Later, we will relax this assumption 
to consider general operations, which can take pure states to mixed states. Such operations 
do involve loss of information. Examples include performing a measurement without retain- 
ing the result, performing an unknown unitary chosen at random, or entangling the system 
with an ancilla which is subsequently discarded. 

The requirement that an operation be local means that the operators U or Mj must have 
a tensor-product structure tJ = tJ ® I , Mj = Mj ® I, where they act as the identity on 
all except one of the subsystems. The ability to use classical communication implies that 
the choice of later local operations can depend arbitrarily on the outcomes of all earlier 
measurements. One can think of an LOCC operation as consisting of a series of "rounds." 
In each round, a single local operation is performed by one of the local parties; if it is a 
measurement, the outcome is communicated to all parties, who then agree on the next local 
operation. 
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B. Entanglement monotones 



For the purposes of this paper, we define an entanglement monotone to be a real- valued 
function of the state with the following properties: if we start with the system in a state p 
and perform a local operation which leaves the system in one of the states pi, • • • , p„ with 
probabilities pi, . . . then the value of the function must not increase on average: 

/(p) > Y.P^f{p,). (3a) 

3 

Furthermore, we can start with a state selected randomly from an ensemble {pkiPk}- If 
we dismiss the information about which particular state we are given (which can be done 
locally), the function of the resultant state must not exceed the average of the function we 
would have if we keep this information: 

Vkpk ■ (3b) 

k \ k J 

Some functions may obey a stronger form of monotonicity, in which the function cannot 
increase for any outcome: 

/(p) > /(P,), Vj, (4) 

but this is not the most common situation. Some monotones may be defined only for pure 
states, or may only be monotonic for pure states. In the latter case, monotonicity is defined 
as non-increase on average under local operations that do not involve information loss. 



C. Infinitesimal operations 

We call an operation infinitesimal if all outcomes result in only very small changes to the 
state. That is, if after an operation the system can be left in states pi, • • • , Pn, we must have 

||p-p,-|| <1, Vj. (5) 

For a unitary, this means that 

il = exp(i£) fti / + i£, (6) 

where £ is a Hermitian operator with small norm, ||£|| <^ 1, £ = . For a generalized 
measurement, every measurement operator Mj can be written as 

Mj = q,{i + e,), (7) 

where < < 1 and Sj is an operator with small norm \ \ej\ \ ^ 1. 

Such measurements are called weak. The term weak measurement, however, is often taken 
to include measurements in which some of the outcomes change the state a great deal, but 
only with very low probability. We do not include such measurements in what follows. All 
outcomes must leave the state almost unchanged. 
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III. LOCAL OPERATIONS FROM INFINITESIMAL LOCAL OPERATIONS 



In this section we show how any local operation that preserves pure states can be per- 
formed as a sequence of infinitesimal local operations. The operations that preserve pure 
states are unitary transformations and generalized measurements. 



A. Unitary transformations 

Every local unitary operator has the representation 

U = e'^, (8) 
where if is a local hermitian operator. We can write 

U= lim(J + zif/n)", (9) 

and define 

i = H/n (10) 

for a suitably large value of n. Thus, in the limit n — oo, any local unitary operation 
can be thought of as an infinite sequence of infinitesimal local unitary operations driven by 
operators of the form 

u,^i + 16, (11) 

where e is a small (||e|| ^ 1) local hermitian operator. 



B. Generalized measurements 



Recently it has been shown 2| that any local measurement can be generated by a 



sequence of weak local measurements. Since a measurement with any number of outcomes 
can be implemented as a sequence of two-outcome measurements, it suffices to show this for 
generalized measurements with two outcomes. 

If the initial state of the system has a density matrix p, the two possible outcomes 
of a measurement with operators Mi and M2 are MipM\/pi and M2pM\/p2-, where 
P12 = Tr(Mi 2pM|2) are the corresponding probabilities. Using the polar decomposition. 



the two measurement operators can be written as Mi = f/i y mIMi and M2 = ^2 y M\M2, 

where JJi and are unitary. Since we have already seen that we can do any local unitary 
transformation by a sequence of infinitesimal steps, if we can first measure the positive op- 
erators sJmIMi and \J MIM2 by a series of infinitesimal steps, we can then apply Ui or U2 

(conditional on the outcome), and can therefore measure Mi and M2 by infinitesimal steps 
as well. So without loss of generality, we consider only positive measurement operators: 
Mj = Mj, i > Mj> 0. Note that in this case. Mi and M2 commute: M1M2 = M2Mi. 

We now decompose this measurement into a series of weak measurements. We can think 
of the procedure as a random walk along a curve in state space; the position on this curve 
is indicated by a single parameter x, with x = being the initial state. The current state 
of the system at any point during the procedure can be written 

M{x)pM{x)/Tt{M\x)p), (12) 
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where 



M(.) = \r^'^^f^-^'^\ .^R. (13) 



In the hmit x — > ±00 the effective operator reduces to Mi,2- 

In it has been shown that depending on the current value of the parameter x, one 
can perform a two-outcome measurement on the system with positive operators M{x, ±e) 
which satisfy 

[M{x,e)]^ + [M{x,-e)]^ = I, 

M(x,±e)M(x) oc M{x±e). (14) 

Since the state is normalized, the factor of proportionality is irrelevant; the two possible 
outcomes simply change the parameter x by +e or — e. Thus the measurement procedure 
is a random walk along the curve M{x), with a step size |e|. We continue this walk until 
|x| > X for some X which is sufficiently large that M{—X) ^ Mi and M{X) ^ M2, to 
whatever precision we desire. It has been shown that the probabilities of the outcomes for 
this procedure are exactly the same as those for a single generalized measurement. The 
exact form of the operators M(x, ±e) is derived in jl^ : 




Mix, ±e) ^ ,IC,^ ,tanh(x ±e)(M|-Mf) 

^ ^ ^' J + tanh(x) (M|-M2) ^ ^ 

where the weights C± are chosen to ensure that these operators form a generalized measure- 
ment: 

C± = {1± tanh(e) tanh(x))/2. (16) 

For any finite x, M{x)~^ is well defined, so from ()14j) and ()15p it is easy to see that if 
|e| <C 1, we have M(x,e) = ^/T/2{I + 0(e)): the measurements are weak. Thus every 
measurement can be implemented as a sequence of weak measurements. Moreover, if the 
original measurement is local, the weak measurements are also local. 

Clearly, the fact that infinitesimal local operations are part of the set of LO means that 
an entanglement monotone must be a monotone under infinitesimal local operations. The 
result discussed in this section implies that if a function is a monotone under infinitesimal 
local unitaries and generalized measurements, it is a monotone under all local unitaries and 
generalized measurements (the operations that do not involve information loss and preserve 
pure states). Based on this result, in the next section we derive necessary and sufficient 
conditions for a function to be an entanglement monotone. 



IV. DIFFERENTIAL CONDITIONS FOR ENTANGLEMENT MONOTONES 

Let us now consider the change in the state under an infinitesimal local operation. With- 
out loss of generality, we assume that the operation is performed on Alice's subsystem. In 
this case, it is convenient to write the density matrix of the system as 

P = PijlmliAliUl ® \ jBC...){'mBC...\, (17) 

i,j,l,m 
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where the set {I^a)} and the set {\jBC...)} are arbitrary orthonormal bases for subsystem A 
and the rest of the system, respectively. Any function of the state f{p) can be thought of 
as a function of the coefficients in the above decomposition: 

f{p)-f{PiJlm). (18) 

A. Local unitary invariance 

Unitary operations are invertible, and therefore the monotonicity condition reduces to an 
invariance condition for LU transformations. Under local unitary operations on subsystem 
A the components of p transform as follows: 

Pijlm — *■ y^'UikPkjpmUjp, (19) 

k,p 

where Uik are the components of the local unitary operator in the basis {|iA)}- We consider 
infinitesimal local unitary operations: 

Uik-{e''),„ (20) 
where £ is a local hermitian operator acting on subsystem A, and 

< 1- (21) 
Up to first order in e the coefficients Pijim transform as 

Pijlm ^ Pijlm ~l~ ^i^i P]ijlm- (22) 

Requiring LU-invariance of f{p), we obtain that the function must satisfy 

J2T^[i,p]ijim^0. (23) 

i,j,l,m^P^^^'^ 

Analogous equations must be satisfied for arbitrary hermitian operators e acting on the 
other parties' subsystems. In a more compact form, the condition can be written as 

iv{|^[£;p]}=0, (24) 

where i is an arbitrary local hermitian operator. 



B. Non-increase under infinitesimal local measurements 

As mentioned earlier, a measurement with any number of outcomes can be implemented 

as a sequence of measurements with two outcomes, and a general measurement can be done 
as a measurement with positive operators, followed by a unitary conditioned on the outcome; 
therefore, it suffices to impose the monotonicity condition for two-outcome measurements 
with positive measurement operators. Consider local measurements on subsystem A with 
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two measurement outcomes, given by operators + M| = /. Without loss of generality, 
we assume 



M2 = \l{I-e)/2, (25) 

where e is again a small local hermitian operator acting on A (in the previous section we 
saw that any two-outcome measurement with positive operators can be generated by weak 
measurements of this type). Upon measurement, the state undergoes one of two possible 
transformations 

P , (26) 

Pl,2 

with probabilities pi.2 = Tr jikf^ 2p|- Since e is small, we can expand 

Ml = -^(/ + £/2-£V8 -■■■), (27) 

M2 = -^(/-^/2-5V8 -■■■). (28) 

The condition for non-increase on average of the function / under infinitesimal local mea- 
surements is 

pJ{MipMi/pi) + P2f{M2pM2/p2) < f{p). (29) 
Expanding in powers of e up to second order, we obtain 

^Tr { f [[^^ P] ^ ^] } + Tr I ^ (xr (.»p - \ {e, p}) | < 0, (30) 

where {e, p} is the anti-commutator of e and p. The inequality must be satisfied for an 
arbitrary local hermitian operator e. 

So long as (jHUIl is satisfied by a strict inequality, it is obvious that we need not consider 
higher-order terms in e. But what about the case when the condition is satisfied by 
equality? In the Appendix we will show that even in the case of equality, (j3(J|) is still the 
necessary and sufficient condition for monotonicity under local generalized measurements. 
There we also prove the sufficiency of the LU-invariance condition This allows us to 

state the following 

Theorem: A twice-differentiable function f{p) of the density matrix is a monotone 
under local unitary operations and generalized measurements, if and only if it satisfies (|24p 
and PH). 

Unitary operations and generalized measurements are the operations that preserve 
pure states. Other operations (which involve loss of information), such as positive maps, 
would in general cause pure states to evolve into mixed states. A measure of pure-state 
entanglement need not be defined over the entire set of density matrices, but only over 
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pure states. Thus a measure of pure-state entanglement, when expressed as a function of 
the density matrix, may have a significantly simpler form than its generalizations to mixed 
states. For example, the entropy of entanglement for bipartite pure states can be written in 
the well-known form Sa{p) = — Tr(pyi log pa), where pa is the reduced density matrix of one 
of the parties' subsystems. When directly extended over mixed states, this function is not 
well justified, since Sa{p) may have a different value from Sb{p)- Moreover, 5'yi(p) by itself 
is not a mixed-state entanglement monotone, since it may increase under local positive 
maps on subsystem A (these properties of the entropy of entanglement will be discussed 
further in section V). One generalization of the entropy of entanglement to mixed states 
is the entanglement of formation which is defined as the minimum of YliPiSAiPi) over 
all ensembles of bipartite pure states {pi,Pi} realizing the mixed state: p = J2iPiPi- This 
quantity is a mixed-state entanglement monotone. As a function of p, it has a much more 
complicated form than the above expression for the entropy of entanglement. In fact, there 
is no known analytic expression for the entanglement of formation in general. The problem 
of extending pure-state entanglement monotones to mixed states is an important one, since 
every mixed-state entanglement monotone can be thought of as an extension of a pure-state 
entanglement monotone. Note, however, that a pure-state entanglement monotone may 
have many different mixed-state generalizations. The relation between the entanglement of 
formation and the entropy of entanglement presents one way to perform such an extension 
(convex-roof extension). For every pure-state entanglement monotone 'm{p), one can define 
a mixed-state extension M(p) as the minimum of 'YliiPi^iPi) o'^sr all ensembles of pure 
states {pi^Pi} realizing the mixed state: p = '^^PiPi- It is easy to verify that M(p) is an 
entanglement monotone for mixed states. On the set of pure states the function M{p) 
reduces to m{p). As the example with the entropy of entanglement suggests, not every 
form of a pure-state entanglement monotone corresponds to a mixed-state entanglement 
monotone when trivially extended to all states - there are additional conditions that a 
mixed-state entanglement monotone must satisfy. On the basis of the above considerations, 
it makes sense to consider separate sets of differential conditions for pure-state and 
mixed-state entanglement monotones. 

Corollary 1: A twice-differentiable function /(p) of the density matrix is a pure- 
state entanglement monotone, if and only if it satisfies (j24j) and (j3(Jj) for pure p. 

For pure states p = \ip){ip\, the elements of p are pijim = c^ijC^imy where the {aij} 
are the state amplitudes: \ip) = ^aij\iA)\jBC...)- Any function on pure states /(p) = /(|V')) 

is therefore a function of the state amplitudes and their complex conjugates: 



By making the substitution Pijim = c^ijC^em i'^^o and (jSOl), we can (after considerable 
algebra) derive alternative forms of the differential conditions for functions of the state 
vector: 



(31) 




(32) 



E 



{sikakj - {e)aij) {Smtain - (^)amn) + c.c. < 0. 



(33) 



dan da, 



mn 
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Here e is a local hermitian operator acting on subsystem A. Analogous conditions must be 
satisfied for e acting on the other parties' subsystems. 



C. Monotonicity under operations with information loss 

Besides monotonicity under local unitaries and generalized measurements, an entangle- 
ment monotone for mixed states should also satisfy monotonicity under local operations 
which involve loss of information. The most general transformation that involves loss of 
information has the form 

Pk . 

where 

Pfc = Tr 1 5^ i,,,p4 ^. I (35) 
is the probability for outcome k. The operators {A^ j} must satisfy 

E = ^- (36) 

We can see that this includes unitary transformations, generalized measurements, and com- 
pletely positive trace-preserving maps as special cases. 

It occasionally makes sense to consider even more general transformations, where the 
operators need not sum to the identity: 

E4A^-^^"- (37) 

This corresponds to a situation where only certain outcomes are retained, and others are 
discarded; the probabihties add up to less than 1 due to these discarded outcomes. We say 
such a transformation involves postselection. 

With or without postselection, we are concerned with the case where all operations are 
done locally, so that all the operators {Akj} act on a single subsystem. Every such trans- 
formation can be implemented as a sequence of local generalized measurements (possibly 
discarding some of the outcomes) and local completely positive maps. In operator-sum 
representation, a completely positive map can be written 

p^5^M,pMt, (38) 

k 

where 

J2MlMk<i. (39) 

k 

Therefore, in addition to (j^^ and (jHUj) we must impose the condition 

f{p)>f[J2MkpMl]. (40) 
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for all sets of local operators {Mk} satisfying (j39p . 

Suppose the parties are supplied with a state pk taken from an ensemble {pk,Pk}- Dis- 
carding the information of the actual state amounts to the transformation 

{pk,Pk} ^ P =^Pkpk- (41) 



As pointed out in [11[, discarding information should not increase the entanglement of the 
system on average. Therefore, for any ensemble {pk,Pk}, an entanglement monotone on 
mixed states should be convex: 



J2pkf{Pk)>f{J2P>^P>^]- (42) 




Condition (j42|) . together with condition (j3(J|) for monotonicity under local generalized mea- 
surements, implies monotonicity under local completely positive maps: 

/ (e ^"'A^^) s Ef*/ s /(p). («) 

It is easy to see that if this inequality holds without postselection, it must also hold with 
postselection. 

It follows that a function of the density matrix is an entanglement monotone for mixed 
states if and only if it is (1) a convex function on the set of density matrices and (2) a 
monotone under local unitaries and generalized measurements. Fortunately, there are also 
simple differential conditions for convexity. A necessary and sufficient condition for a twice- 
differentiable function of multiple variables to be convex on a convex set is that its Hessian 
matrix be positive on the interior of the convex set (in this case, the set of density matrices). 
Therefore, in addition to (j24j) and (j3(Jj) we add the differential condition 

which must be satisfied at every p on the interior of the set of density matrices for an 
arbitrary traceless hermitian matrix a. 

Corollary 2: A twice-differentiable function f{p) of the density matrix is a mixed- 
state entanglement monotone, if and only if it satisfies (j24j) . (j3Up and (j44|) . 



V. EXAMPLES 

In this section we demonstrate how conditions (j21I), ()30|) and (jl^ can be used to verify if 
a function is an entanglement monotone. We show this for three well known entanglement 
monotones: the norm of the state of the system, the trace of the square of the reduced 
density matrix of any subsystem, and the entropy of entanglement. In the next section we 
will use some of the observations made here to construct a new polynomial entanglement 
monotone for three-qubit pure states. 
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A. Norm of the state 



The most trivial example is the norm or the trace of the density matrix of the system: 

h = Tr{p}. (45) 

Clearly Ji is a monotone under LOCC, since all operations that we consider either preserve 
or decrease the trace. But for the purpose of demonstration, let us verify that Ji satisfies 
the differential conditions. 

The LU-invariance condition (|24j) reads 

Tr{^[e,p]}=Tr{[f,p]} = 0. (46) 

The second equality follows from the cyclic invariance of the trace. 

Since the trace is linear, the second term in condition (j3Up vanishes, and we consider only 
the first term: 



Tr 



{^[[^",P],^]}=Tr{[[£,p],£]} = 0. (47) 

The condition is satisfied with equality, again due to the cyclic invariance of the trace, im- 
plying that the norm remains invariant under local measurements. The convexity condition 
(jiij) is also satisfied by equality. 



B. Local purity 

The second example is the purity of the reduced density matrix: 

/2 = Tr{p^}, (48) 

where pa is the reduced density matrix of subsystem A (which in general need not be a 
one-party subsystem). Note that this is an increasing entanglement monotone for pure 
states — the purity of the local reduced density matrix can only increase under LOCC. 

It has been shown in jl^] that every m-th degree polynomial of the components of the 
density matrix p can be written as an expectation value of an observable O on m copies of 
P- 

/(p) = Tr{6p«-}. (49) 

Here we have 



where the components of C are 



Tr{p^}=Tr{c7p^^}, (50) 



Therefore 



Tr||^[f,p]| = Tr|c'([£,p]®p + p®[£,p])} 



TrA {[e, p]aPa + Pa[£, p\a} 

2TrA{pA[£",p]A}, (52) 
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where by Oa we denote the partial trace of an operator O over all subsystems except A. 
If e does not act on subsystem A, then p]a = and the above expression vanishes. If 
it acts on subsystem A, then [e, p]yi = [i,PA\ and the expression vanishes due to the cyclic 
invariance of the trace. 

Now consider condition PU)). If e does not act on subsystem A, then 



From (|3(J|) we get 



[[f,p],£]A = 0. (53) 



1„ ^dh 1 „ [ dH,, /„ 1 



2i:T{[i:T{ep}p-\{e,p}\ \. (54) 



2 



A 



The inequality follows from the fact that {Ti{ip}p — {l/2){e, p})\ is a positive operator. 
If e acts on A, we can use the fact that for pure states 

Tr{p^}=Tr{p|}, (55) 

where B denotes the subsystem complementary to A. Then we can apply the same argument 
as before for the function Tr{p^}. Therefore I2 does not decrease on average under local 
generalized measurements, and is an entanglement monotone for pure states. 

What about mixed states? For increasing entanglement monotones the convexity condi- 
tion ()44|) becomes a concavity condition — the direction of the inequality is inverted. In the 
case of J2, however, we have 

T^I^-^'l-^-t^Klao. (56) 

i.e., the function is convex. This means that Tr{p^} is not a good measure of entanglement 
for mixed states. Indeed, when extended to mixed states, I2 cannot distinguish between 
entanglement and classical disorder. 

C. Entropy of entanglement 

Finally consider the von Neumann entropy of entanglement: 

Sa = -Tr(p^logpA). (57) 

Expanding around pA = I, get 

Sa = -Tr[(p^ - /) + i(p^ - if - ^{pA - 1)' + ...]. (58) 

The LU- invariance follows from the fact that every term in this expansion satisfies ()24j) . If 
we substitute the n-th term in the condition, we obtain 

TT{[i,p]A{pA-ir-') = 0. (59) 
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This is true either because [e, p]a = when e does not act on A, or because otherwise 
[e, p]a = Pa] and the equation follows from the cyclic invariance of the trace. 

Now to prove that Sa satisfies (I^Uj). we will first assume that exists. Then we can 
formally write 

d , dpA d dpA _i 

— log Pa = ——-—log Pa = -J— Pa ■ 

op op OpA op 

Consider the case when i does not act on A. Substituting 5*^ in (jHUI) . we get 



1 

4^' 



dp 



■[[£,p],£]| + Tr 



d 



+ Tr <; ( — ® 

dp 



-Tr 



log Pa 



9p®2 

dpA _ dpA 
dp dp 



TT{ep}p--{e,p} 



TT{ep}p--{i,p} 



_^dpAdpA 
Pa 



Tr{£p}p--{£,p} 



®2 



dp dp 

-Tta |pa' (Tr{fp}p - ^{i, p}^ ^ (Tr{£p}p - ^{e, p} 



-Tr. 



Pa 



1/2 



Tr{£p}p- -{£,p} 



< 0. 



(60) 



(61) 



If p^^ does not exist, it is only on a subset of measure zero - where one or more of the 
eigenvalues of pA vanish. Therefore, we can always find an arbitrarily close vicinity in 
the parameters describing p^, where p^^ is regular and where (jHUjl is satisfied. Since the 
condition is continuous, it cannot be violated on this special subset. 

If i acts on A, we can use an equivalent definition of the entropy of entanglement: 



Sa — S 



B 



-Tt{pb log Pb}, 



(62) 



and apply the same arguments. Therefore 5*^ is an entanglement monotone for pure states. 
The convexity condition is not satisfied, since 



Tr 



d^ 



'Tr{p-/a\} < 



(63) 



This reflects the fact that the entropy of entanglement, like I2, does not distinguish between 
entanglement and classical randomness. 



VI. A NEW ENTANGLEMENT MONOTONE 

It has been shown 3| that the set of all entanglement monotones for a multipartite 
pure state uniquely determine the orbit of the state under the action of the group of local 
unitary transformations. For three-qubit pure states the orbit is uniquely determined by 5 
independent continuous invariants (not counting the norm) and one discrete invariant [l^ 
[13]. Therefore, for pure states of three qubits there must exist five independent continuous 
entanglement monotones that are functions of the five independent continuous invariants. 
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Any polynomial invariant in the amplitudes of a state 



^ aijk...\iA)\3B)\kc) 

i,j,k... 



is a sum of homogenous polynomials of the form [Ji 

P^r-{\^)) = ■ ■ ■ "*„i„fc„...<i,(„)fc,,„,..., (64) 

where a, r, . . . are permutations of (1,2,. . . ,n), and repeated indices indicate summation. A 
set of five independent polynomial invariants for three-qubit pure states is [l^ 



T(AB)C = 


2(1 


-h 


T(AC)B = 


2(1 


-h 


T(BC)A = 


2(1 


-h 


tabc = 


= 2v 





h = Pe,(12) (65) 

h = P(12),e (66) 

^3 = ^(12),(12) (67) 

h = -P(123),(132) (68) 

-^5 = \'^iljlkiOii2j2k2'^i3j3k3^i4j4k4^hi2^i3il^jlj2^j3j4^kik3^k2k4\ ■ (69) 

In the last expression eij is the antisymmetric tensor in 2 dimensions. The first three 
invariants are the local purities of subsystems C, B and A, I4 is the invariant identified by 
Kempe fis'] and is (up to a factor) the square of the 3-tangle identified by Coffman, Kundu 
and Wootters ^19.]. According to [l5| the four known independent continuous entanglement 
monotones that do not require maximization over a multi-dimensional space are 

(70) 
(71) 
(72) 

(73) 

and any fifth independent entanglement monotone must depend on J4. Numerical evidence 
suggested that the tenth order polynomial ctabc = 3 — (/i + /2 + h)h might be such an 
entanglement monotone. However, no rigorous proof of monotonicity was given. Here, 
we will use conditions (|24j) and (j3Up to construct a different independent entanglement 
monotone, which is of sixth order in the amplitudes of the state and their complex conjugates. 

Observe that in (j64|) the amplitudes have been combined in such a way that subsystem A 
is manifestly traced out. By appropriate rearrangement, one can write the same expression 
in a form where an arbitrary subsystem is manifestly traced out. Therefore, any polynomial 
invariant can be written entirely in terms of the components of Tr^i {p} or Ttb{p}, etc. 
This immediately implies that the LU-invariance condition is satisfied, since if e acts 
on subsystem A, we can consider the expression in terms of Pbc.i which, when substituted 
in (j211), would yield zero because [e,p\Bc... = 0. It also implies that in order to prove 
monotonicity under local measurements we can only consider the second term in (j3(J|) . since 
when e acts on subsystem A, we can again consider the expression for the function only in 
terms of psc... and the first term would vanish according to (j53|) . 

We will aim at constructing a polynomial function of three-qubit pure states p which has 
the same form when expressed in terms of pab, Pac, or pBc, in order to avoid the necessity 
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for separate proofs of monotonicity under measurements on the different subsystems. It has 
been shown in that 

h = 3Tr {pab{pa ^ Pb)} - Tr {p^} - Tr 

= 3Tr {pAciPA ®Pc)}- Tr {p^} - Tr {p{} 

= 3Tr {pBciPB ® Pc)} - Tr {p|} - Tr {pi} . (74) 

For local measurements on subsystem C it is convenient to use the first of the above expres- 
sions for I4. The terms Tr{p^} and Tr{p^} are entanglement monotones by themselves. 
This can be easily seen by plugging them in condition (IHUj): 

1 (dTilplr,] 1 id^Tilplr.] ( 1 



dp '^J' M ' 1 Qpm y y r,r 2 

= + 6Tr \^A,B (Tr{£p}p - i{5,p}^^^| > 0. (75) 

These terms, however, are not independent of the invariants I2 and 13. The term which 
is independent of the other polynomial invariants is Tr {p^b(p^ ® Pb)}- When we plug this 
term into condition (j3(Jj) we obtain an expression which is not manifestly positive or negative. 
Is it possible to construct a function dependent on this term, which similarly to Tr |pj4^} 
would yield a trace of a manifestly positive operator when substituted in (jHUIl ? 

It is easy to see that if the function has the form Tr where the operator X{^pab) 

is a positive operator linearly dependent on pab-, it will be an increasing monotone under 
local measurements on C (for simplicity we assume X(0) = 0): 



1 OTr X\pab) ^^Tr X\pab) \ ( 1 x «2 

Tr ^ h% P], ^] ^ + Tr ^ Tr(ep)p - -{5, p} 



4 1 dp LL .rj, J I I y y r,r ^ 

= + 6Tr|x(p^B)X2((Tr{fp}p-i{£,p})^B)} >0. (76) 

Since we want the function to depend on Tx\^pab{pa ® Pb)}, we choose X{pab) = '^Pab + 
Pa ® /fi + /a ® Pb- This is clearly positive for positive pab- Expanding the trace, we obtain: 

Tr ^X^Pab)] = 12Tr{pAB(pA®PB)} + 12Tr{p^5(/A®Pij)} + 12Tr{pi5(pA®/B)} 

+6Tr {pab{Ia ® Pb?} + 6Tr {pab{pa ® Ib?} + 3Tr {pa ® p|} 
+3Tr {p\ ® pb} + Tr {/a ® p|} + Tr {pi ^ Ib} + 8Tr {p^^} 
= 12Tr {pABiPA ® Pb)} + 12Tr {pi^l^A ® Pb)} + 12Tr {plsipA ® /b)} 
+8Tr {p3 } + 8Tr {p|} + 8Tr {p\^} + 3Tr {p^} + 3Tr {pi} . (77) 

One can show that Tr {p^i^ (/a ® Pb)} = Tr {pbc'(pb ® Pc)} and Tr {p^^(pA ® /b)} = 
Tr {pac{pa ® Pc)}- We also have that Tr {pab} = Tr {pc}. Using this and (fTij) . we obtain 

Tr {x'ipAB)} = 12/4 + 16 (Tr {p^} + Tr {p|} + Tr {p^}) + 3Tr {p^} + 3Tr {p|} . (78) 
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This expression is an increasing monotone under local measurements on C. If we add to it 
3Tt {p\q} = 3Tt{pq}, it becomes invariant under permutations of the subsystems. Since 
Tr {pq} is an increasing entanglement monotone, the whole expression will be a monotone 
under operations on any subsystem. We can define the closely related quantity 

= 69 - Tr {{2pAB + Pa®Ib + Ia(E) Pb)'} - 3Tr {pl^} ■ (79) 

This is a decreasing entanglement monotone that vanishes for product states, which is more 
standard for a measure of entanglement. It depends on the invariant identified by Kempe 
and is therefore independent of the other known monotones for three-qubit pure states. 



VII. CONCLUSIONS 

We have derived differential conditions for a twice-differentiable function on quantum 
states to be an entanglement monotone. There are two such conditions for pure-state 
entanglement monotones — invariance under local unitaries and diminishing under local 
measurements — plus a third condition (overall convexity of the function) for mixed-state 
entanglement monotones. We have shown that these conditions are both necessary and 
sufficient. We then verified that the conditions are satisfied by a number of known entangle- 
ment monotones and we used them to construct a new polynomial entanglement monotone 
for three-qubit pure states. 

It is our hope that this approach to the study of entanglement may circumvent some of 
the difficulties that arise due the mathematically complicated nature of LOCC. It may be 
possible to find new classes of entanglement monotones, for both pure and mixed states, 
and to look for functions with particularly desirable properties (such as additivity). There 
may also be other areas of quantum information theory where it will prove advantageous to 
consider general quantum operations as continuous processes. This seems a very promising 
new direction for research. 
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After this paper was written, we became aware of the paper ^23] by Plenio, which claims 
the existence of non-convex entanglement monotones. We would like to point out that such 
a conclusion arises from a different definition of entanglement monotones - namely, functions 
that obey (jSaj) but not necessarily pH]) . At present, we do not know of any argument against 



imposing ()3b|) . which was originally given in and is briefiy discussed in this paper. Since 



ultimately the requirement of ()3b|) is a matter of definition, just like condition (joaj), we have 
conformed to the definition of entanglement monotones given in 

In this version of the paper we have corrected three minor mistakes that appear in the 
published article. The first one was a missing factor of 2 in condition ()30p . the second one 
was in the differential form of the convexity condition (j44p . and the third one was in the 
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proof that the local purity is an entanglement monotone for pure states. Neither change 
alters the conclusions of the paper. 



Appendix: Proof of sufficiency 

The LU-invariance condition can be written as 

where we define 

F(p,e) = /(e-pe--)-/(p) 



(80) 



(81) 



with i being a local hermitian operator. This condition has to be satisfied for every p and 
every e. By expanding up to first order in e we obtained condition (j211), which is equivalent 

to 

f dF{p,i) 
di 



Tr 



1 



e=0 



0. 



(82) 



This is a linear form of the components of e and the requirement that it vanishes for every 
e implies that 

dF{p,e) 



de 



0. 



(83) 



£=0 



This has to be satisfied for every p. Consider the first derivative of F{p,e) with respect to 
Eij, taken at an arbitrary point Eq. We have 



dF{p,e] 



dsi 



dF{p,io + i) 



de 



(84) 



e=0 



But from the form of F{p, i) one can see that F{p, Eq + i) = F{p', i), where p' = e*^°pe 
Therefore 



dF{p,i) 



dF{p\e) 



dsi 



0, 



(85) 

i.e., the first derivatives of F{p,e) with respect to the components of e vanish identically. 
This means that F{p, i) = F{p, 0) = for every i and condition is sufficient. 

The condition for non-increase on average under local generalized measurements (j^^ can 
be written as 

G{p,i)<0, (86) 

where 

G{p,i) = pJ{MipMi/pi)+P2f{M2pM2/P2) - f{p)- (87) 

The operators Mi and M2 in terms of e are given by (j25j) . and the probabilities pi and p2 
are defined as before. As we have argued in section III, it is sufficient that this condition is 
satisfied for infinitesimal e. By expanding the condition up to second order in e we obtained 
condition which is equivalent to 



Tr 



d'G{p,e) 



e=0 



< 0. 



(88) 
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Clearly, if this condition is satisfied by a strict inequality, it is sufficient, since corrections of 
higher order in e can be made arbitrarily smaller in magnitude by taking e small enough. 
Concerns about the contribution of higher-order corrections may arise only if the second- 
order correction to G{p,e) vanishes in some open vicinity of p and some open vicinity of e 
(we have assumed that the function /(p) is continuous). But the second-order correction is 
a real quadratic form of the components of e and it can vanish in an open vicinity of i, only 
if it vanishes for every i, i.e., if 



dsijdski 



= 0. (89) 

£=0 



We will now show that if ()89|) is satisfied in an open vicinity of p, there exists an open vicinity 
of e = () in which all second derivatives of G{p, e) with respect to e vanish identically. This 
means that all higher-order corrections to G{p,i) vanish in this vicinity and ()86p is satisfied 
with equality. 

Consider the two terms of G{p, i) that depend on i: 

G,{p,i)=pif{M,pM,/p,), (90) 



G2{p,i) =P2f{M2pM2/p2). 

They differ only by the sign of i, i.e. Gi{p,e) = 6*2 (p, —£), and therefore 
d^GiipJ) d^G2ip,e) 1 d^Gip,e) 

e=0 



e=0 



If (I89|) is satisfied in an open vicinity of p, we have 

d^Gi{p,e) d^G2{p.e) 



2 deade 



(91) 



(92) 



e=0 



dsijOski 



e=0 



(93) 



e=0 



in this vicinity. Consider the second derivatives of G'(p, i) with respect to the components 
of £, taken at a point Sq: 



d^G{p,e) 



dsijdeu 



d^Gr{p,e) 



e=eo 



dejide 



+ 



d^G2{p,e) 



d'^Gi{p,eo + i) 



dSijOski 



+ 



dsijdski 
d'^G2{p,eo + i) 



i=0 



decide 



e=eo 



e=0 



(94) 



From the expression for Gi{p^e) one can see that e occurs in Gi{p,i) only in the combi- 



nation -v/ ^p\ In Gi(p, Sq + e) it will appear only in J ^ ^i* pJ ^ '1° ^ . But 



'l-i' 



where 



(95) 
(96) 
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So we can write 



^ - ^0 - ^ jl -io-e ^ 1 1 - e' ,^ jl - e' 



where 



and 



Then one can verify that 

Similarly 

where 



Gi{p,eQ + i) =p'Gi{p',i'). 
G2ip,eo + i)=p"G2ip",n, 

i" = i{i + eo)-\ 



p 



I + e^p^I + io If 



Tr 



I + eQP\l I + eq \ . 



(97) 

(98) 

(99) 

(100) 
(101) 
(102) 

(103) 

(104) 



Note that de'^jdsij and de'l^jdsi^ have no dependence on e. Nor do p' and p". Therefore 
we obtain 



5^G(p,£) 



deijdeki 



,d'G^{p',e') 



p 



£=£0 



decide 



+ p 



„ d^G^ip'^e") 



e=0 



dsijdeki 



e=0 



J2 ?fk^p' d'G,ip',i' 



p,q,r,s 



dsij dski 



de' de' 



+ 



J2 ?fk^p" d'G2ip",i" 



p,q,r,s 



dSij dski 



de" de" 

pq^^ rs 



e'=0 



e"=0 



(105) 



We assumed that ()93p is satisfied in an open vicinity of p. If p' and p" are within this 
vicinity, the above expression will vanish. But from (j98|) and (jl03|) we see that as ||eo|| tends 
to zero, the quantities \\p' — p\\ and \\p" — p\\ also tend to zero. Therefore there exists an 
open vicinity of = 0, such that for every io in this vicinity, the corresponding p' and p" 
will be within the vicinity of p for which ()93p is satisfied and 



d'G{p, e) 



deijdeki 



0. 



(106) 



£=eo 



This means that higher derivatives of G{p,e) with respect to the components of e taken at 
points in this vicinity will vanish, in particular derivatives taken at e = 0. So higher order 
corrections in e to G{p,i) will also vanish. Therefore G{p,i) = in the vicinity of p for 
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which we assumed that is satisfied with equahty, which imphes that condition is 
sufficient. 
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